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Rhythmic neuronal activity is a hallmark of many 
higher-level functions of the nervous system including 
the expression of rhythmic motor patterns and genera- 
tion of behaviors such as cognition, memory, arousal, 
and sleep. Many rhythmic networks underlying these 
behaviors consist of one or more pacemaker neurons 
that are coupled with non-pacemaking neurons, and 
exhibit reliable stereotyped activity that however displays 
diverse levels of activity variability between networks - 
and sometimes those circuits that exhibit high levels of 
activity variability even depend upon their variability to 
accomplish some behavioral goal [1,2]. While it is under- 
stood that synaptic feedback from non-pacemaking neu- 
rons coupled to pacemaker neurons can limit variability 
of rhythmic activity [3,4], it is not yet understood how 
specific neuronal dynamics influence variability of the 
coupled network. Here we develop a theoretical approach 
to predict how network activity variability - as quantified 
by the probability distribution of observed periods - 
changes with altered phase response dynamics of non- 
pacemaking neurons, and demonstrate this theoretical 
approach's predictive efficacy in hybrid networks estab- 
lished in the stomatogastric nervous system (STNS) of 
Cancer borealis. These hybrid networks consist of one 
biological pacemaker neuron (Anterior Burster/Pyloric 
Dilator complex [AB/PD]) and one artificial non-pace- 
maker neuron (designed to mimic the burst timing of the 
Lateral Pyloric [LP] neuron), which allows us to vary the 
phase response dynamics of the latter element. 

Our simple theoretical approach to predict the effects 
of network coupling on activity variability is based on 
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three assumptions. First, within the range of variation of 
AB/PD neuron' intrinsic period, it is assumed that there 
is always at least one stable fixed-point in the dynamics 
of the coupled system. Second, it is assumed that any 
drifting of the AP/PD neuron's intrinsic period is signifi- 
cantly slower than the duration of one intrinsic period, 
so that the coupled system is capable of reaching a 
steady-state condition. Third, it is assumed that as the 
period of AB/PD varies, any changes in this neuron's 
phase resetting curve (PRC) are negligibly small. The 
first two assumptions allow us to directly map an AB/ 
PD neuron's intrinsic period onto the network period in 
the coupled system. The last assumption allows us to 
recalculate the map for different intrinsic periods and 
predict the distribution of periods in the coupled net- 
work based upon the known distribution of intrinsic 
periods of AP/PD neuron. Using previously developed 
analyses [5] and a linear ts-tg curve for the artificial 
neuron we are able to qualitatively predict how slope 
and bias of the artificial neuron's ts-tR curve affect the 
distribution of network periods. 

All experiments are performed using standard intra- 
cellular electrophysiological methods for the pyloric net- 
work of the STNS, and hybrid networks are established 
with the dynamic clamp. 
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